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APPLICATION OF INTRINSIC DIFFERENTIATION TO ORBITAL 


PROBLEMS INVOLVING CURVILINEAR 

COORDINATE SYSTEMS 

By James C . Howard 
Ames Research Center 


SUMMARY 


S’- 


62 ^ 


Absolute or intrinsic differentiation has been shown to be an effective 
tool for solving a certain class of vector problems. In many cases, the solu- 
tion of problems involving the rates of change of vectors may be obtained more 
directly by this method than by conventional methods. A single formula with 
sufficient generality to handle a wide range of problems renders application a 
purely mechanical process, requiring little or no ingenuity on the part of the 
analyst. In order to demonstrate its utility, the method has been used to 
process certain fundamental vectors associated with the orbits of planetary 
bodies or space vehicles. The processing of these vectors leads to Lagrange’s 
planetary equations. It is shown that the method has distinct advantages when 
used for obtaining the rate of change of the argument of perifocus of an 
orbiting body. 


INTRODUCTION 




It is well known to students of tensor analysis that problems involving 
the rates of change of vectors of any variance are conveniently solved by the 
method of covariant and intrinsic differentiation. This method of dealing 
with vectors, and the more general entities called tensors, has found little 
or no application to the problems of ordinary vector analysis probably because 
all such problems are amenable to solution by conventional methods. However, 
when covariant and intrinsic differentiation are divorced from the complexi- 
ties inherent in tensor analysis, it is found that certain problems can be 
solved more directly and with much less ingenuity on the part of the analyst 
than is required by conventional methods. This is a consequence of the geo- 
metrical simplification inherent in the method when curvilinear coordinates 
are used. For example, intrinsic differentiation may be used to determine the 
influence of perturbing forces on the time rates of change of certain funda- 
mental vectors associated with the orbits of planetary bodies or space vehi- 
cles. When this is done, a more direct determination of the rates of change 
of some of the orbital elements is possible. In particular, the intrinsic 
derivative of a vector which lies in the orbital plane, and points in the 
direction of the perifocus, leads more directly to a rate of change of the 
argument of perifocus which includes the influence of changes in the longitude 
of the ascending node. Furthermore, since the orientation of a plane is 
uniquely determined by the normal to its surface, the orientation of an orbi- 
tal plane in space is determined by the angular momentum vector. This fact 
may be used to advantage in finding the time rates of change of the elements 


defining the orientation of an orbit plane in the presence of perturbing 
forces. The intrinsic derivative may be used to obtain the rate of change of 
the angular momentum vector and hence the rates of change of orbital plane 
inclination and nodal longitude. These two examples have been chosen to 
illustrate the utility of the method. Moreover, the method may be applied 
with equal facility to any situation in which vector changes are involved. 


NOMENCLATURE 


A 

A 1 (x) 

Ai(x) 

a 

& 

^(x) 

ai(x) 

B-J(y) 

Bj(y) 

b<i(y) 

bj(y) 

e 

"e 

g ij 

gij 

h 

i 

M 

m 

P 


vector 

contravariant vector components in the x coordinate system 
covariant vector components in the x coordinate system 
the semimajor axis 

unit vector which lies in the orbital plane and points in the 
direction of the perifocus 

system of base vectors reciprocal to aj_(x) 
system of base vectors in the x coordinate system 

contravariant vector components in the y coordinate system 
covariant vector components in the y coordinate system 
system of base vectors reciprocal to bj(y) 
system of base vectors in the y coordinate system 
eccentricity 

vector lying in the orbital plane and pointing in the direction 
of the perifocus 

*i * aj 
a 1 • 

angular momentum vector 
orbital plane inclination 
mass of central body 
mass of space vehicle or planet 
perturbing force vector 


2 



r 


position vector 


T 

t 

V 

x 1 


Uj>k] 



ft 

CO 

to 


coordinate transformation 
time 

velocity vector 

components of the position vector in the x coordinate system 
components of the position vector in the y coordinate system 
Christoffel symbol of the first kind 

Christoffel symbol of the second kind 

constant coefficients 

true anomaly 
Kronecker delta 
dynamical constant, G(M + m) 
longitude of ascending node line 

argument of perifocus as measured from ascending node line 
to + ft 


Superscripts 

i,j,k,Z indices of contravariance 

Subscripts 

i,j,k,Z indices of covariance 


ANALYTICAL CONSIDERATIONS 


Transformation Laws for Scalar Components and Base Vectors 


Scalar components .- When referred to a general curvilinear coordinate 
system, a vector A may be expressed in the following form 

A = Ai-a-i 
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If in some expression a certain index occurs twice, this means that the expres- 
sion is to be summed with respect to that index for all admissible values of 

the index, that is, 

9 9 n 

A^a^ = ^ A^a j_ 
i=i 

where A^ are the tensor components of the vector A, and a^, a system of 
base vectors. In accordance with established notation, tensor components will 
be denoted by superscripts and the corresponding base vectors by subscripts. 

In the literature, these vectors are referred to as contravariant vectors, to 
distinguish them from other vectors which are denoted by subscripts. For the 
problems considered in the present report, the distinction between these vec- 
tors disappears. However, it is necessary to keep the distinction in mind, 
because if general coordinate transformations are contemplated the transforma- 
tion law for the components of a contravariant vector denoted by superscripts, 
differs from that for a vector denoted by subscripts. The latter vectors are 
referred to as covariant vectors. For a coordinate transformation T from 
a coordinate system x to a coordinate system y given by 


i 

y 




( 2 ) 


the law of transformation for the components of a contravariant vector A 1 is 
given by (see appendix A and ref. l) : 


J3 J (y) = A 1 (x) (3) 

dx 1 

where A 1 (x) are the contravariant components in the x coordinate system and 
B J (y) are the components when referred to the y coordinate system. For the 
same transformation of coordinates, other vectors, such as the gradient of a 
scalar point function, obey a different transformation law. These are the 
covariant vectors denoted by subscripts. The appropriate transformation law 
for these vector components is (see appendix A) 


B i(y) = A ± (x) (4) 

oy J 

where A^(x) are the covariant components in the x coordinate frame and 
Bj(y) are the covariant components when referred to the y coordinate frame. 
As the following argument shows, the distinction between these two transforma- 
tion laws vanishes when the transformation T is orthogonal Cartesian. Let 
x 1 be the components of a position vector r when referred to the x coordi- 
nate system which is orthogonal Cartesian. Likewise, let yJ be components 
of the same vector when referred to another orthogonal Cartesian system. The 
transformation of coordinates T is given by 
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(5) 


y 1 = 

J 

where the ctj are constants. The position vector r is invariant with 

respect to coordinate transformations. Hence the square of the vector is also 
invariant. Therefore 

x J x j _ yiyi _ aja^x^x k = B^x^x k 

therefore 

a K - 5 k ^ 

where 5 k is the Kronecker delta, that is, (see ref. 2) 

i 1 for j = k 

8 k = 

0 for j / k 

Equation (6) is the orthogonality condition which may be used to solve equa- 
tion (5) for x J . If both sides of equation ( 5 ) are multiplied by a^, 

1 i j i i 
a j a k x = a kT 


therefore 


Therefore, 


J j k i i 
5 £x j = x = a k y 


J 11 

X = a jy 


From equation ( 5 ), it is seen that 


and from equation ( 7 ) 


= rik 


= 4 


It follows from equations (8) and ( 9 ) that 


dy 1 __ 

chcJ Sy 1 
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As a consequence of equation (lO), the distinction between contravariant and 
covariant vectors disappears, when coordinate transformations are confined to 
orthogonal Cartesians systems. This also explains why there is no preoccupa- 
tion with these vectors in the study of ordinary vector analysis. 

Base vectors .- Subscripts assigned to a system of base vectors aq 
indicate that they are covariant in character, and obey the convariant trans- 
formation law. See equation (4) and appendix A. Therefore, if ai(x) are a 
system of base vectors in the x coordinate system, and bj(y) are the cor- 
responding base vectors in the y coordinate system, then 


bj(y) 



( 11 ) 


In this connection it should be noted that to every # system of base vector aq, 
there exists a reciprocal system of base vectors a 1 with the following 
property 


where 


• aq = Sq = aq . a^ (12) 

i 1 for j = i 

“ . 

0 for j f i 


A superscript is assigned to the reciprocal base vectors to indicate their 
contravariant character, and to emphasize the fact that they obey the contra- 
variant transformation law. (See eq. (3) and appendix A.) Hence, _if aA-(x) 
are the reciprocal base vectors in the x coordinate system and b^(y) are 
the corresponding vectors in the y coordinate system, then 


b4(y) = if a x (x) 
ck 1 


(13) 


In a curvilinear coordinate system the base vectors are, in general, not unit 
vectors, but are functions of the coordinates; that is, 


aq = aq (xU 2 , . . ., x n ) 

(14) 

a4 ss a^x^-x 2 , . . x n ) 

(15) 


The base vectors may be obtained as follows: let dr be the differential of 

a position vector r and let dx 1 be the corresponding differentials of the 
position components. Then by substituting dr for A, and dx^ for A 1 in 
equation (l), we have 


dr = dxAaq 


(16) 
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From equation (l6) the base vectors aq are given by 


a± = 


a? 

ax 1 


( 17 ) 


In an orthogonal Cartesian frame of reference, the base vectors aq consti- 
tute a triad of mutually orthogonal unit vectors, that is, vectors of unit 
length. However, in problem formulation, it is usually convenient to use a 
more general curvilinear coordinate system. When this is done, the magnitudes 
of the base vectors generally differ from unity. 


Vector Derivatives and the Christoffel Symbols 

The scalar product of any two base vectors aq and aj may be defined as 
follows : 


a j “ Sij a j • a i 


(18) 


Likewise, the scalar product of the reciprocal base vectors a 1 and may be 
defined as 


a* 3 = g 1 * 3 = a J 


— 1 
a 


( 19 ) 


The symmetry of gqj and gij follows from the nature of the scalar product. 
Certain combinations of the partial derivatives of the scalar products with 
respect to the system coordinates are useful in obtaining the derivative of a 
vector, or in formulating the equations of motion in a general curvilinear 
coordinate system. The definitions that follow are ascribed to Christoffel 
and are called Christoffel symbols (see ref. 3)- There are two of these 
symbols, the first of which is defined as 


[ij, k] = ! 




( 20 ) 


The Christoffel symbol of the second kind is 




[ij. 


l] 


( 21 ) 


The utility of the Christoffel symbols is immediately apparent when an attempt 
is made to find the partial derivative of a base vector, or its reciprocal, 
with respect to any system coordinate* Any vector A may be expressed in the 
form of equation (l). Furthermore, since the base vectors are in_general 
functions of the coordinates, it follows that the derivative of A with 


7 


respect to any coordinate must involve the Christoff el symbols. From equa- 
tion (l), the partial derivative of the vector A with respect to the 
coordinate x^ is given by 


However, since ap 


Likewise, 


and 


Since 


it follows that 


SA _ SA^ — , A i Sa^ 

3x k = 3x k 1 Sx* 


a j “ Sij 

hA = 


Sgp/j Sap 


dx^ dxk 


aj + ap 


Sgpv Sa,. 


x — 


= — r • a k + ap 


dx<J SxJ 


a. = 


Sr 

dx 1 


da,. 


Saj 

dxk 


= Saj . a k + a i • ^ 
dx^ Sx^- dx^ 

ik 


Sa^ 

dxj 


_S_ f Sr 

dxi VdxJ 


ox j dxj \Sxi 

From equations (23) through (26) 

Sap _ 

— . a k = [ij, k] 
oxJ 


da., 

dxi 


( 22 ) 

(23) 

(24) 

(25) 

(26) 
(27) 


From equation (12), the rate of change of the base 
to Xi assumes the form 


Sap 

s^j 


[ij, k]a k 


vector 


ap with respect 


(28) 


Equation (28) gives the required rate of change of the base vector, with 
respect to a system coordinate, in terms of the Christoffel symbol of the 
first kind and the reciprocal base vectors. A more convenient form is 
obtained if both sides of equation (28) are multiplied scalarly by the recip- 
rocal base vector a^ to yield 

— — . = [ij, k]a k . a 1 (29) 

Sx J 

From equation (19), it is seen that 

_k -l k l 

a • a = g 
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therefore 


(30) 


da ± -1 

7 • a 


[ij f k]g k * 


In terms of the defining formula (2l), equation (30) may be rewritten as 
follows: 



Therefore, 



By substitution of equation (32) in equation (22) the partial derivative of a 
vector A with respect to the system coordinate x^ is 


& 


Mi % 

Sxk 1 


+ A 1 



( 33 ) 


The indices i and l in the second term on the right side of equation ( 33 ) 
are dummy 1 indices, and may therefore be replaced by any other convenient 
indices. In order to have a common base vector aq, equation ( 33 ) may be 
rewritten as follows 

II = fat + AO 

c)xk \3xk 


dA dx k 

5Ai dx^ 

c)x k dt 

the intrinsic derivative, or the derivative with respect to time, may be 
obtained from equation (3*0 in the following form: 


Furthermore, since 


and 



dA _ 
dt 



( 35 ) 


x As already indicated, a repeated index implies summation with respect to 
that index. Since the summation index can be changed at will, it is usually 
referred to as a dummy index. Of course, the range of admissible values of 
the index must be preserved. 
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In an orthogonal Cartesian reference frame 


s iJ 


= a. 




= g 


ij 


Therefore, since all these scalar products are constants, it follows that the 
Christoffel symbols vanish. In this case, the covariant derivative, (3*0 
reduces to the sum of the partial derivatives of the components along a set 
of fixed axes 


5A _ dA^ — . 
dx^ dx^ 


i - 1, 2, 3 


Likewise, the intrinsic derivative of a vector reduces to the ordinary time 
rates of change of the components along a set of fixed axes. 


For a general space of three dimensions, equation (35) assumes the form 



The formidable looking equations ( 36 ) through (39) f or intrinsic deriva- 
tive of a vector in a general space of three dimensions contain 27 
Christoffel symbols. Because of the symmetry of the Christoffel symbols. 


< 


and the number of independent Christoffel symbols reduces to 18. Furthermore, 
for the three-dimensional spaces most commonly used, equation ( 36 ) reduces to 
a manageable form. If a base vector of unit length be denoted by aq, then 
in a cylindrical coordinate system, 

\ 

a-i = ax , gii = 1 

a 2 = x 1 ^ , g 22 = (xl) 2 > (41) 

— A 

a 3 = a 3 j g33 = 1 

> 

As a consequence of equations (4l) , there are only two nonzero Christoffel 
symbols in a cylindrical coordinate system, embedded in a space of three 
dimensions . These are 



Hence, a vector referred to this coordinate system has a time rate of change 
as follows: 



In three-dimensional spherical coordinates, 

ai = a x , gn = 1 

a 2 = x x a 2 , g 22 = (x 1 ) 2 

a 3 - xi sin x 2 a 3 , g 33 = (x 1 sin x 2 ) 2 

In this case there are six nonzero Christoffel symbols. These are 



(44) 


(45) 


11 


When these values of the Christoffel symbols are substituted in equation (35)> 
the time rate of change of a vector, referred to a three-dimensional spherical 
coordinate system, assumes the following form: 



APPLICATIONS 

Intrinsic Derivatives of Certain Fundamental Vectors Associated 
With the Orbits of Planetary Bodies and Space Vehicles 

Certain fundamental vectors associated with an orbit in space, and which 
are subject to change in the presence of perturbing forces, provide a useful 
area of application for intrinsic differentiation. One example of such a 
vector is the angular momentum vector h which lies in the direction of the 
normal to the orbital plane. This vector may be used to determine the orien- 
tation of an orbital plane in space. Furthermore, the rate of change of 
this vector in the presence of perturbing torques may be used to determine 

the rates of change of orbital plane 
inclination and nodal longitude. 

The intrinsic derivative of a vector 
e, which lies in the orbital plane 
and is directed to the perifocus, 
may be used to determine the rates 
of change of the argument of peri- 
focus and the eccentricity. A 
third vector may be used to complete 
the orthogonal triad. 

The angular momentum vector . - 
Equation (35) may be used to obtain 
the rates of change of the angular 
momentum vector, and hence the 
rates of change of the orbital ele- 
ments defining the orientation of 
an orbital plane in space. In 
spherical coordinates, equation (35) 
assumes the form shown in equa- 
tion (46). The coordinate uses are 
chosen as follows (see sketch (a)): 



Sketch (a) 
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1. _The ai axis is taken to be coincident with the angular momentum 
vector h. 

2. The a 2 axis is in the direction of increasing polar angle, that 
is, in the direction of increasing orbital plane inclination. 

3. The third axis £3 is in the direction of the ascending node line 
and completes the mutually orthogonal triad of axes. 


In this coordinate system the base vectors are given by equations (4k) . Let 
the vector A appearing in equation ( 46 ) have components as follows 


A 1 = h 


Ap — Ap = 0 


(47) 


When substitutions are made from equations (44), (45), and (47) in equa- 
tion (46), it is found that 


dA 

dt 


dh 

dt 


dh 

dt 




1 _ 

.x 1 


dx 2 

dt 


i A 

x a 2 


_1 

xi 


therefore 


dh _ f dh 
dt \dt 


a x + h 


dx 2 a 

dT a2 


+ h sin x 2 


dx~ 

dt 


dx 3 

dt 


sin x 


,A I 

a3 


& 


(48) 


These are the component rates of change of the vector h, referred to the x 
coordinate system with base vectors aj_. However, it may be more convenient 
to use some other reference frame. Assuming that the perturbing forces are 
referred to the radial and transverse direc- 
tions in the plane of the orbit, and in the 
direction of the normal to the orbital 
plane, the vector dh/dt should be trans- 
formed to a reference frame having these 
directions. If the coordinates in this ref- 
erence frame be denoted by y 1 , where i= 1 , 

2 , 3, a coordinate transformation Ti must be 
determined before the components of dh/dt 
can be referred to the y coordinate frame. 

Let 
the 


bj(y) be a system of base vectors in 
y coordinate system (see sketch (b)). 


In terms of the notation already estab- 
lished, the transformation of base vectors 
assumes the following form 


£1 = -sin (7 + w)a 2 + cos(y + co )& 3 
b 2 = -cos(y + to)S 2 - sin(y + u )& 3 
b 3 = t x 


(49) 



Sketch (b) 
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If A i (x) are the components of dh/dt in the x coordinate system and 
BJ(y) are the corresponding components in the y coordinate system, the 
transformation law, equation (3) gives 

B'ky) = hi A^x) 

Sx 1 


However, 


fiJ -q 

OX 


See appendix A. Therefore, equation (3) may be written in the following 
alternative form 


B J (y) = (b i • ai)A 1 (x) 


From equation (48) 


A 1 = — , A 2 = h , A 3 = h sin x 2 (! 

Substituting from equations (51) in equation (50), and using equations (49) 
we obtain the transformed components of the vector dh/dt as follows 

2 3 \ 

B 1 (y) = -h sin(y + to) ^4- + h C os (7 + to) sin x 2 


B 2 (y) = -h cos (7 + to) 


h sin(7 + to) sin x 2 

CLu 


However, 


b3 ^> - 1 


^!=fxP 


where r is the position vector of the planet_or space vehicle in the y 
coordinate system with base vectors bj, and P is the perturbing force 
vector. Therefore, 


r X P = rbi X (bjbj) * P 


r X P = rb x X (bib! + b 2 b 2 + b 3 b 3 ) * P 


r X P = r(b 3 b 2 - b 2 b 3 ) * P 
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From equations (52) and (53)> it follows that 


h cos(7 + to) sin x 2 — h sin(7 + to) 

-.3 2 

h sin (7 + to) sin x 2 -r— + h cos (7 + to) —— 

dt w dt 

dh 


6yf 


0 

r(b 3 


These equations give 


1 ■ *■>« • 


dx 2 _ r cos(7 + to) ^ . p) 


dt 


h 


dx 3 _ r sin(7 + co) (£ 3 . p) 
^ h sin x 2 


P) 


( 5*0 

(55) 


Equation ( 5*0 gives the rate of change of inclination, and equation (55) the 
rate of change of nodal longitude. 

The orbital vector e .- Con- 
sider the vector e, which lies in 
the orbital plane and is directed 
to the perifocus. (Bee sketch (c) 
and ref. 4.) In a cylindrical ref- 
erence frame with coordinates 
denoted by the base vectors 

aj_ have the following values: 


ax = &i 

a 2 = x x a 2 

&3 — 


(56) 


In this case, it should be noted 
that x 2 is the nodal longitude. 

In the case previously considered, 
x 2 was the orbital plane inclina- 
tion. In this reference frame, 
which is chosen because the influ- 
ence of nodal longitude and orbital plane inclination appear in the formula- 
tion, the vector e has components as follows (see appendix B and sketch (c)) 



that is. 


e = e l cos wai + s - 1 - n - - -- os — a 2 + sin co sin ia.3 


e = A^a-j_ = A x ai + A 2 a 2 + A^a^ 
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where 


A 1 = e cos co 


A 2 = (e sin co cos i)/x x l 
A 3 = e sin co sin i 


(57) 


In cylindrical coordinates, equation ( 35 ) assumes the form shown in equa- 
tion ( 43 ). Substituting from equations ( 42 ) and ( 57 ) in equation (43) gives 


de __ f 

dt I U 


(e cos to) - xi^f • e sin co cos 1 
dt dt x 1 


^1 


d f e sin to cos i \ . 1 f 

_dt \ x 3 - J xi \ 


+ dt s ^ n 10 s ^ n i )^ 3 


e cos w + e ^ ■ si . n - aj . c . os 1 
dt dt xi 


x 1 ^ 


Therefore, 

de 
dt 


dx* 


e cos co - e sin toco - e —— sin to cos i } ai 

dt 1 


+ I e sin co cos i + eco cos co cos i - e — sin to sin i + e cos to ) a 2 
' dt dt ' 


+ ( e sin to sin i + eto cos to sin i + e — sin to cos i ) a 3 

V dt / 3 J 


( 58 ) 


These are the component rates of change of the vector e referred to the 
coordinate system with base vectors aq. As in the case of the angular 
momentum vector, it may be more convenient to use some other reference frame. 
Assuming again that the perturbing forces are referred to the radial and 
transverse directions in the plane of the orbit, and in the direction of the 
normal to the orbital plane, the vector de/dt must be transformed to a ref- 
erence frame having these directions. A coordinate transformation T 2 must 
be determined before the components of de/dt can be referred to this coor- 
dinate frame. Let the base vectors in this reference frame be again denoted 
by bj, then the base vectors bj and Sq are related as follows: 


bx = cos(y + w)£x + sin(y + co) cos ia 2 + sin(y + co) sin ia 3 
b 2 = -sin(/ + <o)ax + cos(y + co) cos ia 2 + cos(y + co) sin ia 3 1 

A 

b3 = -sin ia 2 + cos ia3 


(59) 


Corresponding to the transformation of coordinates T 2 , the transformation 
law for the components of the vector de/dt is given by equation (50) . In 
this case, the components Ai(x) appearing in the right side of equation (50) 
may be obtained from equation (58) . They are 
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A 1 = 


e cos co - ew sin to 


&xr 

dt 


sin to cos 1 


•\ 


A 2 s* ( e sin co cos i + ew cos co cos i - e — sin co sin i + e cos co 
\ dt dt 

A 3 = sin co sin i + eco cos co sin i + e ^ sin co cos i^ 


di - 


dx 2 


> ( 60 ) 


When substitutions are made from equations (59) and (60) in equation ( 50 ), it 
is found that the vector de/dt has the following components in the direc- 
tions of the base vector bj 


B 1 (y) = ( e cos 7 + eel) sin 7 + e sin 7 cos i 


dt 


B 2 (y) = ^-e sin 7 + eto cos 7 


+ e 


dx 2 

dt 


cos 7 cos i 


di 

dt 


dxf 


B (y) = ( e — sin co - e sin i cos co 


J 


( 61 ) 


In order to express the rate of change of the vector e as a function of the 
perturbing force vector P, it is necessary to examine the equation of motion 
of a particle in a noncentral, gravitational force field. It is shown in 
equation (B19) that 


de 

dt 


2h 


L P 


— bib 2 - ( — sin 7b 2 b 2 + — b 2 bi ) - 


\ 


re 


/ 


sin 7b 3 b 3 


( 62 ) 


By equating components from equations (6l) and (62), and remembering that in 
this case dx 2 /dt = dft/dt, we obtain the following equations: 


e cos 7 + eco sin 7 = ~ (b 2 * P) - e 4^ sin 7 cos i ( 63 ) 

eco cos 7 - e sin 7 = - sin 7b 2 + ^ bj^ • P - e 4j^ cos 7 cos i (64) 

e 4^~ sin co - e ^ sin i cos co = - — sin 7(b 3 * P) ( 65 ) 


The rate of change of the eccentricity e and the rate of the change of the 
argument of perifocus may be obtained by solving equations ( 63 ) and (64) . 
When this is done, it is found that 


“ = - ■— cos 7(bi • P) 


+ — P- sin 7(2 
eh 


+ e cos 7)(b £ 


73 \ cL£2 

p) - it cos 


( 66 ) 
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It is seen that the influence of changes in the longitude of the ascending 
node is given by the third term on the right side of equation (66), that is. 


bh 


-cos i 


(67) 


However, if the argument of perifocus as measured from the inert ially fixed 
X 1 axis is denoted by co, then 


CO rs (jj + ft 


(68) 


In this case, the total contribution from the rate of change of the nodal 
longitude is 


From equation ( 55 ) 



ft(l - cos i) 


2ft sin 2 



ft 


(b 3 • p> 

h sin 1 


(69) 


therefore 


• _ r 

bh h 


sin(y + co)tan i^(b 3 • P) 


( 70 ) 


Substituting from equation (70) in equation (68) gives 


A* 

CO 



bi 


+ sin 7(2 + e cos 
eh 



sin(y + co)t 



P 


( 71 ) 


The derivation of this result should be compared with the approach used in 
reference 4 , where spherical geometry had to be used to determine the influ- 
ence of inclination and nodal longitude. In that case, the problem had to be 
solved in two parts, whereas in the present case the problem is solved in one 
step without appeal to system geometry, subsequent to the choice of curvi- 
linear coordinates. 


Likewise, on solving equations (63) and ( 64 ) for e, it is found that 


® ~ "{( 77 sin 7 ) b i + 


k[ cos 


7(2 + 


e cos 


7) + 




(72) 


If required, the rate of change of the semimajor axis may be obtained with the 
aid of equation (72) and the following relationship 


h 2 = qa(l - e 2 ) 
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therefore 

therefore 


therefore 


2h — = p(l - e 2 )a - 2paee 
at 


a = ■■ 2 ^ e . 


e + 


2h dh 


(i(l - e 2 ) |i(l - e 2 ) dt 


„ _ 2ae • j. 2hr r* . ^ 

a = — e + -7- -T- (b 2 * P) 


(l - e 2 ) n(l - e 2 ) 


Therefore , 


a = 


2af 

h 


e sin 7bi + (l + e cos 7)b 2 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Feb. k, 1965 


APPENDIX A 


TRANSFORMATION FORMULAS FOR THE BASE VECTORS AND THEIR RECIPROCALS 


The transformation formulas and, hence, the covariant or contravariant 
character of the base vectors and their reciprocals may be obta in ed as 
follows: Let the differential of a position vector be denoted by_ dr. Then 

if ai(x) are the base vectors in the x coordinate system, and bj(y) the 
base vectors in the y coordinate system, the differential dr may be 
expressed in the following alternative forms 


dr = aj^xjdx 1 = bj(y)dy^ = bj(y) dx 1 


dx 1 


therefore 


Likewise, 


therefore 


a^x) = dd. bi(y) 

dx 1 J 


af(x) 1^- dy J = bj(y)dy J 


(Al) 

(A2) 




b i(y) = a ± (x) 

3 W 


(A3) 


It is seen from equations (A2) and (A3) that the base vectors ap and bj 
obey the covariant transformation law; consequently, the use of subscripts 
is justified. 


Reciprocal Base Vectors 


To each system of base vectors a^ there exists a reciprocal system of 
vectors aJ with the following property 


^i * a J = 8^ = a j • aq 
where 5 ^ is the Kronecker delta, that is, 

1 for j = i 


(A 4 ) 


S'? = 


0 


for j ^ i 


Scalar multiplication of each side of equation (A2) by b^(y) gives on using 
(A 4 -) (see ref. 5) 


b J (y) • ai(x) = 

dx 1 


(A5) 
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Similarly, from equations (A3) and (A 4 ) it is seen t ha t 


a 1 (x) • bj(y) = — — 7 

V 


(A6) 


Equation (Al) referred to the reciprocal system of base vectors assumes the 
form 

dn = a i (x)dx i = bJ(y)dy^j 


therefore 


therefore 


and 


therefore 


dyj = bj(y) • a 1 (x)dx i 
^i = b J (y) • a jL (x)dy 1 


dxi = dyj 


From equations (A7) and (A8) 


dx 1 


dxi 


therefore 


F-Cxjax^ = bJ(y) dx ± 
oy d 


al (x) = b j(y) 

dyJ 


Likewise, from equations (A7) and (A9) 


bJ(y) = H a i (x) 

ox 1 


Cat) 


(A8) 


(A9) 


(A 10 ) 


(All) 


From equations (AlO) and (All), it is seen that the reciprocal base vectors 
a 1 , obey the contravariant law of transformation; therefore, the superscript 
notation is justified. 


Vector Transformations 

Equations_(A 10 ) and (All) may be used to obtain the transformation law 
for a vector A, where 


A = A-^a^ = Aja^ 


(A 12 ) 


21 



If A = Ai(x)aj_(x) when the vector A is referred to the x coord ina te 
system, and if _A = BJ(y)bj(y) when referred to the y coordinate system, 
invariance of A requires that 

B J (y)bj(y) = A 1 (x)a i (x) (A13) 

From equations (A2) and (AI 3 ), the appropriate transformation law is obtained 
as follows: 


B J(y) = |4 Ai(x) (A14) 

OX 1 

Equation (Al4) is the contravariant transformation law for the components of 
the vector A. When A is referred to the x coordinate system with base 
vectors aj_(x) , which obey the covariant transformation law, the components 
Ai(x) obey the contravariant transformation law, and, hence, the use of 
superscripts is justified. If A is referred to the reciprocal base system 
a 1 , then from equation (A12) 


A = A-^a 1 

On a transformation of coordinates from the x coordinate system to the y 
coordinate system, the invariance of A requires that 

A i (x)a 1 (x) = Bj(y)b^(y) (A 15 ) 

From equations (A10) and (A15), the appropriate transformation law is obtained 
as follows 

Bj(y) = 14 A i (x) (Al6) 

oyJ 

It is seen that when a vector A is referred to a coordinate system by use 
of reciprocal base vectors, which obey the contravariant law, the correspond- 
ing components of A obey the covariant law, and the use of subscripts is 
therefore justified. 
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APPENDIX B 


EQUATION OF MOTION OF A PARTICLE IN A NONCENTRAL GRAVITATIONAL FORCE FIELD 


The equation of motion of a particle of unit mass moving in an inverse- 
square law, central force field is 


a 2 ? _ a * 

at 2 " " r 2 bl 


(Bl) 


Vector multiplication of each side of equation (Bl) by the angular momentum 
vector h gives 


a r w r 


at c 


x h = - JL (bi x h) 


h = r X V — (i^yjba 

On substitution from equation (B3) in equation (B2) it is seen that 


(B2) 

(B3) 


therefore 


x h = (ti7)b 2 

at 2 


a /v x h) _ a * 
at VT7 "at 1 


(B4) 


The integral of equation (B4) is given by 


V X h * . - 

= bi + e 


(B5) 


where ¥ is a constant vector of integration. The vector "e may be 
expressed in terms of its scalar magnitude and a vector of unit length as 
follows : 


e = ea 


where §, is a vector of unit length. See sketch (c). Substituting equa- 
tion (B 6) in equation (B5) gives 


V X h 

a 


= bi + ea 


(b6) 


(bt) 


Equation (B7) may be solved to obtain the position vector r. Scalar multi- 
plication of each side of equation (B7) by r gives the following equation 
for r: 


V X h 


= r + e(r * §.) 
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Therefore 


and 


( r x y . ) — — = r(l + e cos 7) 
V- 

r - th. 

1 + e cos 7 


(B8) 


It is seen that e, the magnitude of the constant vector of integration 
appearing in equation (b 6 ), is the orbital eccentricity. Vector multiplication 
of each side of equation (B7) by h gives 


Therefore, 


h X (V X h) = |jh X (bi + ei.) 




+ e(b3 X a) 


(B9) 


With the notation of sketch (c), the unit vector & may be expressed in the 
following form: 

A A 

a = (cos 7)bi - (sin 7)b 2 (BIO) 


If the assumption of an inverse -square law central force field is not satis- 
fied, the equation of motion must be modified accordingly. In the presence 
of a perturbing force P, the equation of motion becomes 


dV p _ JJ£ 

dt r 3 


(Bll) 


Furthermore, in the presence of the perturbing force vector P, the assumption 
of constancy no longer applies to the vector e. Hence, 


de _ JdV 
dt “ ]dt 


X h + V X (r X P) - ^ [(r X V) X r] 


} 


(B 12 ) 


Therefore, 


and 


de f &V L |ir 


^ dt V dt + r 3 


X h + V X (r X P) 


|i - = P X h + V X (r X P) 

dt 


(B 13 ) 


The first term on the right side of equation (B 13 ) may be written in the fol- 
lowing alternative form: _ 

P x h = h(bib 2 - b 2 bi) • P (Bl 4 ) 


2 k 


Likewise, 


r X P = r(b 3 b 2 - b 2 b 3 ) • P 


(B15) 


When substitutions are made from equations (B9) and (B15), the second term on 
the right side of equation (B13) assumes the form 

V X (r X P) = ip -jbib 2 + e|ab 2 + (b 2 • a)b 3 b 3 Jj- • P (Bl6) 

Substituting for a from equation (BIO) in (Bl6) gives 

V X (r X P) - [(l + e cos 7)bib 2 - e sin 7(b 2 b 2 + b 3 b 3 )] • P (B17) 

h 

From equations (Bl4) and (B17) it follows that 

PXh+Vx(rXP)= j^(2h)bib 2 - sin 7b 2 b 2 + hb 2 bi^ - pp sin 7b 3 b 3 • P 

(B18) 

Therefore, 


de (2h * * [re 

at = ij 12 " L~ S1 


sin 7(b 2 b 2 + b 3 b 3 ) + ^ b 2 bi | [• • P 


*]} 


(B19) 
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